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Abstract 

The photon vacuum polarization effect in curved spacetime leads to birefringence, 
i.e. the photon velocity becomes c±5c depending on its polarization. This velocity 
shift then results in modified photon trajectories. We find that photon trajectories 
are shifted by equal and opposite amounts for the two photon polarizations, as 
expected by the sum rule p]. Therefore, the critical circular orbit at u = 1/3M in 
Schwarzschild spacetime, is split depending on polarization asu= 1/3M ± AS(M) 
(to first order in A), where A is a constant found to be ~ 1CT 32 for a solar mass 
blackhole. Then using general quantum modified trajectory equations we find 
that photons projected into the blackhole for a critical impact parameter tend 
to the critical orbit associated with that polarization. We then use an impact 
parameter that is lower than the critical one. In this case the photons tend to 
the event horizon in coordinate time, and according to the affme parameter the 
photons fall into the singularity. This means even with the quantum corrections the 
event horizon behaves in the classic way, as expected from the horizon theorem 0. 
We also construct a quantum modified Schwarzschild metric, which encompasses 
the quantum polarization corrections. This is then used to derive the photons 
general quantum modified equations of motion, as before. When this modified 
metric is used with wave vectors for radially projected photons we obtain the 
classic equations of motion, as expected, as radial velocities are not modified by 
the quantum polarization correction. 
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Chapter 1 
Introduction 

1.1 General Relativity 

Since the birth of special relativity, in 1905, the nature of space and time has 
been demoted to a relative entity, known as spacetime, which is stretched and 
contracted depending on an observer's frame of reference, while the speed of 
light, c, has taken the pedestal of an absolute and universal speed limit, 
unaffected by any transformation of reference frame. From this emerged a 
generalized theory of relativity, which portrayed the gravitational field in a 
new and revolutionary way: where it didn't depend on a propagating field 
but on the nature of spacetime itself. In this view matter (or energy) is 
said to curve and modify the surrounding spacetime, this then results in 
photons and particles tracing out shortest paths between two points, known 
as geodesies. Therefore, gravitational forces become a manifestation of the 
curved spacetime due to the presence of matter [HE]. In this general rela- 
tivistic framework spacetime is described by the metric and the motion 
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of particles are described by the interval equation: 

> Time-like (c < 1) 
k 2 = g^Wk" =0 Light-like (c = 0) (1.1) 

< Space-like (c > 1) 

where /c M = ^r- 1 . For flat spacetime or a local inertial frame (LIF), where 
g^u is replaced by the diagonal Minkowski metric rj^ = (1,-1,— 1,-1), the 
interval equation becomes: 

,9 , „,„ dt 2 dx 2 dy 2 dz 2 ,„ . 

for c = 1. Apart from resolving the problems associated with Newtonian me- 
chanics, such as describing the perihelion advance of Mercury, the strongest 
aspect of general relativity was its predictive power. One of its most radical 
claims, and the building blocks of the theory itself, was that gravitational 
fields affect radiation, which was then confirmed through the observation 
of starlight deflection by the sun. From this emerged some profound and 
fantastic possibilities such as black holes and gravitational lensing. The 
gravitational effect on light rays also leads to the possibility that photons 
could follow stable orbits around stars (discussed in chapter EJ), and it's this 
possibility in which we will be interested. 

1.2 QED in a Curved Spacetime 

Even though photon trajectories are modified in a curved spacetime, and the 
resulting curved paths are described by general relativity, this bending of 
light was, for a long time, considered to have no effect on the velocity of the 

photon. This view shifted slightly when, in 1980, Drummond and Hathrell[T] 

1 For photons this becomes fe M = written in terms of the affine parameter A 
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Figure 1.1: First order, a, vacuum polarization Feynman diagram contribut- 
ing to the photon propagator 

proposed that a photon propagating in a curved spacetime may, depending on 
its direction and polarization, travel with a velocity that exceeds the normal 
speed of light c. This change in velocity would then result in trajectories 
other then the ones described by "classical" general relativity. This effect is 
simply described as a modification of the light cone in a LIF: 



Where a is the fine structure constant and a /JiU (R) is a modification to the 
metric that depends on the Riemann curvature at the origin of the LIF. This 
correction is seen to arise from photon vacuum polarizations in a curved 
space time, Fig. II . H Qualitatively it can be thought of as a photon splitting 
into a virtual e + e~ pair, so at the quantum level it is characterized by the 
Compton wavelength A c ; then, when this quantum cloud of size 0(X C ) passes 
through a curved spacetime its motion would be affected differently to that 
described by general relativity, possibly in a polarization-dependent way 

This effect of vacuum polarization is considered through the effective 
action: 



k 2 = T)» v Wk v = 



(rj^ + aa^(R))k^k 



(1.3) 



S — S\ + 5*2 



(1.4) 
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where Si is the Maxwell electromagnetic action in curved space time: 

Si = ~J d A x^F^F^ 

= d^A u -d u A^ (1.5) 

and S 2 is the part of the action that incorporates the effects of virtual electron 
loops in a background gravitational field. As we are only concerned with 
the propagation of individual photons it must be quadratic in A^, and the 
constraint of gauge invariance then implies that it must depend on F^ v rather 
than Ap. Also, as the virtual loops give the photon a size of A c , S 2 can be 
expanded in powers of A 2 = m~ 2 , thus the lowest term in the expansion 
would be of order m -2 , which is the term corresponding to one electron loop. 
With these constraints there are only four independent gauge invariant terms, 
which can be chosen to be: 

S 2 = J d^V^giaRF^F^ + bR^F^FZ 

+ cR^F^F^ + dD^ v D a F°) (1.6) 

The first three terms represent a direct coupling of the electromagnetic field 
to the curvature, and they vanish in flat-spacetime. The fourth, however, 
is also applicable in the case of flat-spacetime, and represents off-mass-shell 
effects in the vacuum polarization. In the values for a, b, c, and d have 
been determined to 0(e 2 ). The constant d is obtained by comparing the 
coefficient of the renormalized flat-spacetime photon propagator associated 

with the Feynman diagram in Fig. II . if to the result of the same order given 

2 The photon propogator with the vacuum polarization in flat-spacetime is given (in the 
Feynman gauge) by: -> + jclpu, where 1^ = (rj^q 2 - q^q u ){l - ^2^2 + • • •) 
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by the effective action S; and a, b, and c are obtained by comparing the 
coefficients of the coupling of a graviton to two photons 3 to the same result 
obtained from S 2 . In this way the constants are given as: 

5 a 26 a 2 a 24 a 

a — 6 = c= a = (1.7) 

720 7T 720 7r 720 tt 720 tt v ' 

Then, as the equations of motion for the electromagnetic field are given by: 

SS (1.8) 



using the modified action (|1.4j) we find: 

D,F^ + ^ = (1.9) 

From this we can see that D^F^ is of 0(e 2 ), therefore, the term with coef- 
ficient d in Eqn. (|1.6jl will be of 0(e 4 ), hence we can omit it from the final 
equation of motion. In this way we find 0: 

- ^{2bR^D»F x » + Acg^R^D^F^} = (1.10) 

which is the Maxwell equation in curved spacetime, incorporating the cou- 
pling of curvature with vacuum polarization effects. Using this modified 
Maxwell equation and the methods of geometric optics (described in Chap- 
ter EJ) it is possible to derive the quantum modified light cone and geodesic 
equations. Then, using these, we are able to determine the quantum modified 

trajectories in curved spacetime. 

3 Deduced from the matrix element < 7(92, l3)T^ v ^{qi, a) >, where is the energy 
momentum tensor, pQ 
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1.3 The Equivalence Principle and Causality 

The equivalence principle exists in two forms: weak and strong. The weak 
equivalence principle states that at each point in spacetime there exists a local 
Minkowski frame, which is a fundamental requirement of general relativity 4 . 
The strong equivalence principle (SEP), on the other hand, states that the 
laws of physics are the same in all LIFs at different points in spacetime, 
and at the origin of each LIF they take the special relativistic form. Then 
the coupling of curvature to the electromagnetic field in the effective action, 
Eqn. (|1.4j) . is a violation of the SEP. Due too this violation of the SEP, QED 
in curved spacetime remains a causal theory - despite the modification to the 
physical light cone. 

This can be seen more clearly by considering Global Lorentz invariance 5 
(GLI), which is the special relativistic equivalent of the SEP. In special rela- 
tivity GLI states that faster than light signals automatically imply the pos- 
sibility of unacceptable closed-time-loops, Fig. 11.21 That is, if you can send 
a signal backwards in time in one frame, then it should be possible in any 
frame. However, if you break this GLI, a signal backwards in time in one 
frame does not automatically imply you can send a signal backwards in time 
from any frame. Therefore, in the case of QED in curved spacetime, due to 
the breakdown of the SEP we retain the fundamental property of causality: 
faster than light signals can be seen to go backwards in time in a certain 

frame, but this no longer implies that you can send a signal backwards in 

4 This implies that general relativity is formulated on a Riemannian manifold 
5 Global Lorentz invariance states that the laws of physics are the same in all inertial 
frames. 
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Figure 1.2: An unacceptable closed-time-loop. 

time from any frame. Thus, in this case, spacelike motion does not necessarily 
imply a causality violation 6 . 

1.4 Critical Stable Orbits 

One of the simplest curved spacetimes is the (Ricci-flat) Schwarzschild space- 
time, which describes the spherically symmetric geometry outside a star. The 
geometry of such a spacetime structure is given by the line interval: 

ds 2 = (1 - ^)dt 2 - (1 - —)- x dr 2 - r 2 d6 2 - (r 2 sin 2 6)d<p 2 (1.11) 

In this spacetime the equations of motion, derived for k 2 = 0, have a special 

solution for null geodesies. This solution describes a light ray in a circular 
6 Further discussion of causality is given in ^H] 
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orbit with a radius r = 3M and an associated critical impact parameter. 
Therefore, a photon projected from r — > oo with the critical parameter tends 
to the stable circular orbit by spiralling around it. But, in the context of 
vacuum polarization effects, we may derive the null equation of motion us- 
ing the quantum modified Maxwell equation. In this way the orbit equation 
should give us stable circular orbits that are dependent on the photon polar- 
ization, i.e. splitting the circular orbit depending on polarization. Therefore, 
a photon coming in from r — > oo, with a particular polarization and impact 
parameter, should tend to its associated critical stable orbit by spiralling 
around it in the classic way. 



Chapter 2 

Null Dynamics in Schwarzschild 
Spacetime 

In this chapter we will give a brief overview of null dynamics in Schwarzschild 
spacetime. Starting with the geodesic equations and the line interval for pho- 
tons we will derive the orbit equation, which will then be used to determine 
the critical stable orbit for photons and the associated impact parameter. 
We will then go on to show that decreasing the impact parameter, from the 
critical value, causes the photon trajectory, given as a function of 0, to spiral 
into the singularity, but as a function of t it tends to the event horizon. The 
main aim of this chapter is to familiarise ourselves with the classical solutions 
of photon trajectories in Schwarzschild spacetime, so we can then compare 
the equivalent results for the quantum modified case. 

2.1 Equations of Motion 

When we consider motion in a plane, where 9 = | and 9 = 9 = 0, the 
geodesic equations for Schwarzschild spacetime, (jA.10j) - (|A.13j) . can written 
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as: 



2M . ■■ 2M . . . 

^-—) t+ — H = (2.1) 



M . 9 , 2iW i M . 2M X , 9 . . 

— t 2 - r0 2 +1 f - t(1 ^ = 2.2 



2rr<\> + r 2 = (2.3) 

And the spacetime line interval becomes: 

2M, .„ , 2M, . 9 9 . , 

(1 - — )t 2 - (1 - — )^r 2 - r 2 2 = X (2.4) 

where .fT = for photons and ±1 for space-like or time-like motion. Now, 
to derive the orbit equation ^ we use Eqns. (|2.1j) . (|2.3j) and the interval 
equation (12.4(1 . in the simplified form: 

= i + — (2.5) 
2r • 

= 0+— (2.6) 

r 

2ilf 9 , 2M n9 . 9 9/ 2M - 9 , . 

K(l = -r 2 + 1 2 t 2 -r 2 l U 2 2.7 

In order to specify photon trajectories we can set K = at anytime, however, 
then r would be interpreted as an affine parameter and not proper time. To 
proceed we divide (Q by f - and $M b Y f-' 



i" 2Af, 2ilf , 

• + — 1 v 

t r z r 



= I + )-V (2.8) 



These can written as: 



2r , , 

= + — 2.9 

fli r 



= — (lnt + ln(l - — )) (2.10) 
r/ • 2r 

= -(ln0 + -) (2.11) 

dr r 
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Then, solving these we have: 

i = (2.12) 

r 

= J - (2.13) 

Where E and J are constants of integration denoting total energy and an- 
gular momentum about an axis normal to the plane 6 = tt/2. Now, using 
Eqns. (I2.12|) and (|2.13|) to substitute for t and <j) in Eqn (|2.7jl and rearranging, 
we have: 



r 2 + r 2 1 — + Kl = 1 2 1 X E 2 2.14 



Rearranging and simplifying, this becomes: 



2 



(-) 2 + (l 2.15 

dr r r z 

,dr. ^ r , 2M,,K J 2 ..i , 
(_) = £[! -(!-_)(_ + _)], (2.16) 

We now have the components of the four momentum P Q in spherical polar 
coordinates 



P a = (t,r,M) 

= [1-^ + ^)1^0,^) (2.17) 

where D = 4 is the impact parameter, F = (1 — — ) and K = for photons 
and F = +1 for particles. 

2.1.1 Orbit Equation 

In order to derive orbit equations that are physically understandable we need 
to represent them as r(<f>), r(t) and 4>(t). In this form we can analyse the 
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orbit paths as a function of rotational angle and the time taken to reach a 
certain point along the angle. Therefore, by using 



dr dr dcf) dr J 
dr deb dr dcbr 2 



(2.18) 



we can rewrite (I2.15B as 



dr I 2 9M T 2 

e 2 = Qr- + (i-—)(K + J -) 

d(p r r z 

(_) 2 = (E 2 -K)— + 2MK— -r 2 + 2Mr (2.19) 

Now, transforming u = r" 1 , so at r = oo we have u — 0, Eqn. (|2.19|) becomes: 

AV = ( E 2 -K)^ + 2MK^-r 2 + 2Mr 
dm J A J z 



,du. (E 2 — K) u 

{-^f = - — j 2 — - + 2MK— — u + 2Mu (2.20) 



Doing similar manipulation for (f>(t) we have 

d(p d<p dr 
dt dr dt 



(2.21) 



and using Eqns. (|2.12|) and (|2.13j) in this, and transforming u = -, we have: 

dd> D, 2M 9/ . 

f t = -(1 - — ) = Du\l - 2Mu) (2.22) 

Finally, we can determine u{t) by inserting: 

du du dt du, „ 9 , , , m i 

— = — = — L>u 2 (l - 2Mu " 2.23 

d<b dt deb dt [ v ;J v 7 



into Eqn. (|2.2flj) . which gives 



A 2 = [Du 2 (l - 2Mu)] 2 [ ^ E \ K) + 2MK^- — u 2 + 2M« 3 ] (2.24) 
dt J 2 J 1 
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Now, the Eqns. (127201) . (127221) and (1272311 are the general equations that de- 
termine photon and particle trajectories in the plane 6 = ix /2. If we take 
K = we then have the required photon trajectory equations 

A 2 = -L - u 2 + 2M« 3 = /(«) (2.25) 



' 2 (l-2Mu)] 2 [ 1 
[Du 2 (l - 2Mu)] 2 f(u) (2.26) 



;^) 2 = [Du 2 (l-2Mu)} 2 [—-u 2 + 2Mu 3 



We can also write another equation, specifically for null radial geodesies. By 
using the fact that = 0, Eqn. 1)2.131) implies J = 0, then Eqns ()2.12jl and 
1)2.15)1 become 

2M s dt 

1 )^ = E 2.27 

^- = ±E (2.28) 
ar 



Combining these we have: 



dr , 2M . 

S = ±(l-— ) (2-29) 



2.2 Orbits in Schwarzschild Spacetime 



We will now discuss the solutions of Eqns. (12725ft . (12361) and fj2729ft . We will 
start off with the simple case of the radial geodesic and then go onto the case 
of the general orbits. For the general case we will solve 1)2.2 51) and 1)2.26)) 
to determine the critical stable orbits and the associated impact parameters. 
We will then go on to study the trajectories of photons as they fall into the 
singularity 
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r 



Figure 2.1: Radial Null geodesies take an infinite coordinate time to reach 
the event horizon at 2M 



2.2.1 Null Radial Geodesic 

Eqn. (|2.29j) can be solved by rewriting it as: 

dt , 2AI . , 

Tr = * l --r l ^ (2-30) 

t = ±(r + 2M\og(^—-l)) + constant± (2.31) 

This solution (specifically the — one) gives the trajectory of a photon coming 
from infinity and into the black hole. It shows that the photon takes an 
infinite coordinate time to reach the horizon, which can be seen in Fig. 12.11 
However, if we solve Eqn. (l2~28l) . i.e. in terms of the affine parameter, we 
can show that the photon reaches and crosses the event horizon without ever 
noticing it, 

r = ±Et + constant± (2.32) 
this can be seen in Fig. 12.21 
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2.2.2 General Null Geodesies and Critical Orbits 



We will first solve Eqn. (I2.25jl in order to determine the the critical stable 
photon orbits. In order to do this we first consider the point of equilibrium 
and the associated impact parameter. This equilibrium point occurs when 

du 

H = (2 - 33) 

which means as the photon is orbiting the black hole u does not change; so 
if the radial distance does not change it implies a circular orbit. Therefore 
we must solve 

f(u) = ^-u 2 + 2Mu 3 = (2.34) 
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The sum and product of the roots u\, u 2l and u 3 of this equation are given 
by 1 

1 1 

ui + u 2 + u 3 = — and u lU2 u 3 = ~ 2M£)2 (2-35) 

This shows that f(u) = must have a real negative root, and the two 
remaining roots can be real (distinct or coincident) or be a complex conjugate 
pair; however, the occurrence of coincident positive real roots implies the 
existence of a circular orbit. Thus, if a coincident root occurs it should be at 
the point given by the derivative of f(u) 

f(u) = 6Mu 2 -2u = (2.36) 

Which then has the solution U\ = u 2 = (3M) _1 . For this solution the impact 
parameter of equation ()2.34|) is D = (3\^3)M. From the product condition, 
equation ()2.35|) . we find that the roots of f{u) = are 

u 1 =u 2 = — - and u 3 = -— — and D = (3\/3)M (2.37) 
6M blVl 

Therefore, when the impact parameter is D = (3v^3)Af then % vanishes for 
u = (3M) _1 , which implies a circular orbit of radius 3M is an allowed null 
geodesic [2]. 

Now we can consider a photon at u = with an impact parameter D = 
(3a/3)M. This, then, gives a trajectory of a photon spiralling in and tending 
to the critical orbit at u = (3M) _1 . The general differential equation for 
this impact parameter is given by rearranging and substituting for D in 
Eqn. (|2~55j) 

{ p = mu + JL )iv -JLr (2.38) 



1 for au 3 + bu 2 + c = we have u\ + 1*2 + 1*3 = — 77 and U1M2M3 
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From |2j we have the solution to this as 

where 0o is a constant of integration, given by: 

tanh 2 (-^>o) = ^ (2.40) 

which gives: u = (r — ► oo) when = 0, and u = when <fi — > oo. 
Therefore a null geodesic arriving from infinity with an impact parameter 
D = (3v3)M approaches the circle of radius 3M, asymptotically, by spi- 
ralling around it, as can be seen in Fig. 12. 3f . Also, numerically solving 
Eqn l|2.2fij) we can show that as time increases u tends to which can be 
seen in Fig. 12.41 

Finally we can show that when an impact parameter other than D = 
(3\/3)M is used, for example if we set D = (3y/3)M — 0.1, the solution 
of Eqn. (|2.25j) shows that the photon falls past the critical orbit, through 
the event horizon, and into the singularity, Fig. 12.51 Using this new impact 
parameter in Eqn. (|2.26|) and, again solving numerically, we see that the 
photon comes in from infinity and tends to the event horizon, u = ^g, 
asymptotically in time t, Fig. 12.61 



2 This figure was plotted using Eqn Q2.39I1 . We also obtained the same plot by numeri- 
cally solving Eqn Ij2.38l) in Mathematica. 
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Figure 2.4: For impact parameter D = (3v^3)M a null trajectory tends to 
u = ox? asymptotically with time. 
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Figure 2.6: Null geodesic, with impact perimeter D = (3\/3)M — 0.1, arrives 
from infinity and tends to u = Jrj asymptotically in time t. 



Chapter 3 

Quantum Gravitational Optics 



3.1 Photon Propagation in Curved Spacetime 

Maxwell's equations, in curved spacetime, 

= D^ u (3.1) 

= D^F uX + D u F^ + D x F^ u (3.2) 
=S> = d^A v - d v A„ (3.3) 

cannot be solved explicitly, and even in cases of extreme symmetry explicit 
solutions are difficult; this is due to the fact that as curved space acts as 
a dispersive material (i.e. bending light rays) plane wave solutions do not 
exist. 

3.1.1 Geometric Optics in Curved Spacetime 

In dispersive materials, where light rays are bent, we can consider the solution 
of Maxwell's equations to be a simple perturbation of the plane wave solution. 
For example in curved space, relative to an observer, the electromagnetic 
waves can appear to be plane and monochromatic on a scale that is much 

20 
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larger compared to the typical wavelength, but very small compared with 
the typical radius of curvature of space time. Such "locally plane" waves can 
be represented, in geometric optics, by approximate solutions of Maxwell's 
equations of the form[6j: 

T» v = Re(Fr + teFf + . . .)e^ (3.4) 

and the electromagnetic field vector, defined by Eqn. (13. 3 J) , takes the form: 

A" = Re(A^ + ieA% + . . .)e*f (3.5) 

where the electromagnetic field is written as a slowly- varying amplitude and a 
rapidly-varying phase. The parameter e is introduced in order to keep track 
of the relative order of magnitude of terms, so in curved space Maxwell's 
equations can be solved order-by-order in e. In this formulation the wave 
vector is then defined as the gradient of the phase of the field, = D^O = 
d^d, which in terms of the quantum interpretation is identified as the photon 
momentum. We can also write A^ = Aa 11 , where A represents the amplitude 
and (normalized —1) specifies the wave polarization. These 

vectors then satisfy the condition = 0. 

Geometric Optics and Null Dynamics 

In this notation Eqn. 1)3.1)) can be written, to leading order as: 

<9,(FfVf) = Fre^C-d.O) 
=> K F x = (3.6) 
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and Eqn. (|3.3|) becomes: 

F Xilv e^ = (dpAiv - d v A ly )e^ 

= l -(d,9A lu - d u 9A lfi y 8 - 
F lflv = k„A lu - k v A lfl (3.7) 

Now, by combining these, we have: 

k,Fr = k tl (k»A1-k»A l t) = 
= (k^k^A" — k u (k^A^) 
=> k 2 a u = (3.8) 

and from this we can deduce that k 2 = 0, i.e. k^ 1 is a null vector. Also, it 
follows from the definition of W as a gradient that D^ky = D u k^, so 

FD^r = WD"^ = \® v k 2 = (3.9) 

Using this, and the fact that light rays are defined as the curves given by 



s) where ^ = k^ we can derive the geodesic equation as follows 



= WDyfr" 

d 2 x v dx^ dx x 

= -jf + ^irir (3 - 10) 

3.2 Quantum Modified Null Dynamics 

As was seen in Sec. II . 2| using the effective action, Eqn. (|1.4j) . the equation 
of motion becomes: 

= - ^{2bR^D»F x » + Acg^R^D^F^} (3.11) 

and the Bianchi identity, Eqn. (|3.2j) . remains unchanged. Now, as before we 
can determine the quantum modified light cone and geodesic equations. 
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Quantum Modified Light Cone 

Substituting Eqn. flJQJ in ED) we fin d, again to O(i): 



+ 4^- J R AtrA ^( J D^)e J fF 1 A "] 

= fc„Fr - [2bR^F^ + ^^fe"^ (3.12) 
and using Eqn. (|3.7jl we have: 

+ Acg UT R pT \ p k p {k x A{ - k p A x )} 

using k^A^ = this becomes: 

= k^Al - -LfibR^^Al-k'A*) 

+ Acg VT R liTXp k^(2k x A p )} (3.13) 

where A\ = A\a v , and the last line is simplified by relabeling of indices. 
Now, Contracting with Aa v and eliminating A we have: 

= knk^ H -Rn\k^ k -Rn T \ p k^k &^a!l 
m z f ml 

This then gives the quantum modified light cone: 

k 2 - ^R pX k p k x - ^R prXp k p k x a T a p = k 2 - 6k(a) = (3.14) 

where we have replaced the constant c = —3^: with c a = for convenience 
of interpretation, i.e. the sign of the light cone is immediately obvious from 
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R pTXp k p 'k x a T a p . As we will be working in the Schwarzschild spacetime the 
quantum modified light cone for the Ricci flat case {R^x = 0) is given by: 



k 2 = -^R prXp k»k x a T a p = Sk{a) (3.15) 



Here the sign of the light cone depends on polarization and the photon tra- 
jectory; if the correction is positive we have space-like motion, and if it's 
negative we have time-like motion. 

Quantum Modified Geodesic Equation 

The photon trajectories corresponding to the quantum modified equation of 
motion, 1)3. llj) . can be represented by a generalised version of Eqn. 1)3 



= ^[k 2 - ^-R pX k p k x -^R prXp k p k x a T a p 
2 m l mi 



= -D u k 2 - —D u [bR pX k p k x + Ac a R pTXp k p k x a T a p ] 



2 " m 2 
d 2 x v „„ dx p dx x 



= ^r + T 



ds 2 mA ds ds 

dx^ dx^ 

- —d v [{bRfr + ±c a Rp„ 1T a a a T )——] (3.16) 

where we have used /c M = ^f-, and covariant derivative in the second term is 
replaced by a partial derivative as it's acting on a scalar. In Ricci spacetime 
this equation becomes: 

3.3 Horizon Theorem and Polarization Rule 



There are two general features associated with quantum modified photon 
propagation^. First, it is a general result that the velocity of radially di- 
rected photons remains equal to c at the event horizon. Second, for Ricci 
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flat spacetimes (such as Schwarzschild[Tj and KerrjH] spacetimes), the veloc- 
ity shifts for the two transverse polarizations are always equal and opposite. 
However, this is no longer true for non-Ricci flat cases (such as Robertson- 
Walker spacetime[T]). In these cases, the polarization averaged velocity shift 
is proportional to the matter energy-momentum tensor. These features can 
be easily shown by using the Newman-Penrose formalism: this characterises 
spacetimes using a set of complex scalars, which are found by contracting 
the Weyl tensor with elements of a null tetrad[2j. 

Newman-Penrose Formalism 

We choose the basis vectors of the null tetrad as[Zj: Z M = the photon mo- 
mentum. Then, we denote the two spacelike, normalized, transverse polar- 
ization vectors by a M and \f and construct the null vectors m M = ^(a M + z6 M ) 
and m M = 775 (a^ — ih^). We complete the tetrad with a further null vector 
n^, which is orthogonal to m fl and rh fl . We then have the conditions: 

I - m — I ■ m — n - m — n ■ fh — (3.18) 
from orthogonality, and: 

I ■ I — n ■ n — m ■ m — m- fh — (3.19) 
since the basis vectors are null. Finally, we impose: 

I ■ n = —m -m = l (3.20) 
The Weyl tensor, given in terms of the Riemann and Ricci tensors, is: 



1 



C'fivyS R^ivyS r. {jl^RuS VffRfiS VfiS-Rvy ^]v&R 



2 

g(%y?/V<5 - V^sVu^R (3.21) 
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where R^ = rj vS R^ ul s and R = rj^R^; and the Weyl tensor satisfies the 
trace-free condiation: 

V^C^s = (3.22) 

and the cyclicity property: 

Cl234 + C*1342 + C*1423 = (3.23) 

Now, using the null tetrad, we can denote the ten independent components 
of the Weyl tensor by the five complex Newman-Penrose scalars: 



$ = 


-C^s^mTm 6 


*1 = 


-C^n u Vm s 


#2 = 


-C^ l5 l»m u m<n s 


^3 = 




\P 4 = 


—C^sn tJ, fh u n y rh s 



(3.24) 

3.3.1 Polarization Sum Rule 
Ricci Flat Spacetime 

In Ricci flat spacetime, summing the quantum correction over the two po- 
larizations leads to the following polarization sum rule: 

^<5&(a) = (3.25) 

a 

This can be proven by suming the quantum correction in Eqn. (|3.15j) over 
the two polarizations, 

Sk{a) = R^sk^^a'a 5 (3.26) 
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In the Newman-Penrose basis, using k = I, a = ^(m+fh), b = — ^-{m— m), 
and the fact that C^ ul s = R^s for the Ricci fiat case, we have: 

^R^skWa'a 6 = ^C^sFPim" + m v )(m 5 + m s ) 

a 

- ^C^sPP (m u - fh v ) (m s - m 5 ) 

= C^ s PP{m u rh 5 + rh u m 5 ) (3.27) 

This particular contraction is equal to zero as it's not part of the complex 
scalars in Eqns. (|3.24|) : hence the sum of the two quantum corrections is zero. 
This implies the trajectory (and velocity) shifts are equal and opposite. 

Non-Ricci Flat Spacetime 

For the non-Ricci flat spacetimes the polarization sum rule is given as: 

6k(a) = — %(2b - 8c a )Tn,k ft k v (3.28) 

a e 

where T^ v is the energy-momentum tensor. This can be shown by proceeding 
as before, but now we include the Ricci tensor and scalar, as in Eqn. (|3.21|) : 

R^sk^k"<a u a s = C^ 1& PP{m u m 5 + m v m s ) - R^PP (3.29) 

As before the first term on the RHS is zero, and the second term is only 
dependent on the photon momentum. Then, combining this with the Ricci 
term in Eqn. (13. 14 J) we have: 

5k(a) = — (26 - 8c a )i?^Fr (3.30) 

a e 

Finally replacing the Ricci tensor with the energy-momentum tensor, by 
using the Einstein equation 

R»v = -8nT^ u + -Rg^v (3.31) 
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we obtain Eqn. (l3~28l) . 

3.3.2 Horizon Theorem 

At the event horizon, photons with momentum directed normal to the hori- 
zon have velocity equal to c, i.e. the light cone remains k 2 = 0, independent 
of polarization^. 

This can be easily proven for the Ricci flat spacetime, using the orthonormal 
vectors k a = (E k ,E k ,0,0), a b = (0,0,1,0) and a b = (0,0,0,1). Therefore, 
using these vectors in Eqn. (|3.15jl . we have: 

k 2 = ^R abcd k a k c a b a d = (3.32) 

So, in Ricci flat spacetime all radially projected photon trajectories remain 
unchanged. 

It is also possible to prove the horizon theorem for the general case (for 
Ricci and non-Ricci flat spacetimes) that the light cone at the event horizon 
is unchanged. This can be seen in the null tetrad, so that the physical, space- 
like, polarization vectors a M and 6 M lie parallel to the event horizon 2-surface, 
while k^ is the null vector normal to the surface. Then, from Eqn. ()3.14j) we 
have for the two polarizations: 

k = —Rn-yk^k^ -\- — Rni/^gk^ k^ CL CL 




Using Eqn. (13.311) and the fact that C [m v m & + m u m 5 ) = 0, we can 
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write this as: 

k 2 = -^(26 - 4c a )T„JT ±^C m /r(mV + m v m 8 ) (3.34) 
and in terms of the Newman-Penrose scalars this can be written as: 

k 2 = -^(2&-4 Ca )r,/r±^*o, (3.35) 

where the simplification in the last term on the RHS is possible as is 
real for Schwarzschild spacetime. In general Eqn. (|3.35jl is non zero, how- 
ever, at the event horizon the terms: T^IH 1 and \l/o are zero for stationary 
spacetimes0 [H] 1 . 

Physically the Ricci term represents the flow of matter across the horizon 
and the Weyl term represents the flow of gravitational radiation[2j, and both 
are zero in classic general relativity; and as, even with the quantum modi- 
fication, the light cone remains unchanged at the event horizon, means the 
event horizon is fixed and light cannot escape from inside the black hole. 



1 Stationary spacetimes are independent of time and may or may not be symmetric 
under time reversal 



Chapter 4 

Quantum Modified Trajectories 



In this chapter we will analyse how the classical null trajectories, described in 
Chapter El are modified in Schwarzschild spacetime, due to quantum modifi- 
cations of the equations of motion of general relativity. Using Eqn. (|3.14j) we 
will calculate the quantum corrected version of Eqn. 1)2. 25|) . which will then 
describe the quantum modified motion of a null trajectory in a Schwarzschild 
spacetime. Using this, and following the methods of Chapter O we will begin 
by studying simple critical circular orbits at the stationary point, u = 1/3M. 
As stated in the polarization rule, the critical orbit should be shifted up and 
down by equal amounts, depending on the polarization of each photon. Also, 
we will show that these modifications are only valid if the "classic" impact 
parameter corresponding to the stationary orbit is adjusted, depending on 
the photon's polarization, in order to compensate for the orbit shift. This 
information, of the modified orbits and the corresponding impact parameter 
for each polarization, will then be used to determine the general trajectory 
of a (vertically or horizontally polarized) photon coming in from infinity and 
tending to one of the two shifted critical orbits. We will then go on to show 
that these quantum modifications have no effect on the event horizon, that 

30 
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is, when the impact parameter is accordingly adjusted and a photon falls into 
the singularity the event horizon remains fixed at u = 1/2M. 



4.1 Quantum Modified Circular Orbits 

In Schwarzschild spacetime, using k v = dx u /dr, Eqn. ()3.15|) can be written 
as: 

= r 2 - (1 - — )V + r 2 (l - — )6 2 

+ (l-—y^ + ^(l-—)R^ s kWa s (4.1) 



and Eqn. (ITTTI) as: 



r ' m 2 , ' r 



o = t+— (i-— r 1 ^ 

^(fl^aV^S (4.2) 
m~ ar ar 



+ ^-^L^AR^Ja^-f) (4.3) 
r z r mf ar ar 



2r0 4c Q _ , „ ,. r fix' 4 aV , 
= + — --ffy {R lunS a v a s — — 4.4 
r mf ar ar 



4c„ „ , ,, aoV* ax 7 , 



= -^^(^Wa^^r^) 



e 



where c a = 3^-. Now, in order to consider quantum modified circular orbits 
we require three things: (i) the Riemann tensor components, (ii) the photon 
wave vectors, k^ = for circular orbits, and finally (iii) the photon polar- 
ization vectors, a M . Due to the circular nature of the orbit the simplest basis 
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to work in is the orthonormal basis. In this frame the required polarization 
and wave vectors, for circular orbits, are simply given as: 

a M = (0, 1, 0, 0) Planar Polarized (4.6) 



a M = (0, 0, 1, 0) Vertically Polarized (4.7) 
k^ = — — = (E k , 0, 0, E k ) 4-momentum along (4.8) 

GST 

Now using these photon vectors, the quantum modifications in Eqns. (|4.1|) - 
(|4.5j) . can be expanded to give: 

• For the planar polarized case we have: 

Rabcdk a k c a a = R arcr k k = E^{R aTCT k fc°) 

= E\[R tr tr + R<t>rcf>r + ^Rtr<f>r) (4-9) 

• For the vertically polarized case we have: 

R ab cdk a k c a b a d = R adce k a k c = E 2 k (R aece k a k c ) 

— El(R m9 + R^e^e) (4-10) 

Using the six independent components of the Riemann tensor in the orthonor- 
mal basis 1 , we find: 

2M 2M 

Rtrrt — —Rtrtr — -jj - =^ Rtrtr — V*-H) 



1 Which we have calculated in Appendix A 
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From Eqns. 1)3.15)1 . ()4.16|1 and ()4.18jl we can see that: 



3E 2 M 

k 2 ~ — \ — For Planar Polarization (4.20) 

3E 2 M 

k ~ \ — For Vertical Polarization (4.21) 



This implies that, as the light cone for planar polarization is positive, it 
represents a photon trajectory with a speed < c, and as the light cone for 
vertical polarization is negative, it represents a photon trajectory with a 
speed > c. 

Planar Polarization 

Considering the planar polarized case first, we can use Eqns. (|4.16j) and 1)4.17]) 
to rewrite Eqns. P~Tf - ((4~4]) as: 

9 , 2M N 9 . 9 9/ 2M ■„ 
= r 2 -(l ) 2 t 2 + r 2 (l )0 2 

2Af 8c a/ 2M X 3E 2 M 

+ 1 r 2 2 + ^l 4.22 

r m~ r r 6 

.. 2M , 2M . . . 
= t + — (1-— )~Vt (4.23) 



M . 2M . 9 , 2AT 
= r - — (1 - — ) _1 r 2 - r(l 



+ — 1 * 2 + Hf^- 4.24 



2r<p 

+ — (4.25) 



where Eqn. 1)4.5)1 becomes zero. Now, using the solutions of (]4.23)1 and 1)4 .2 5)1 . 
as given by 1)2.121) and 1)2.13]) . we can rewrite Eqn. 1)4.22]) as a simple quantum 
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modified trajectory equation for circular orbits with radius r and in a plane 



9 = -- 2 

u 2 • 



, 2 ^ 2 „ 2M.J 2 8c QM 2M.3E 2 M 
r-E 2 + (l )— H r(l 



To make this equation more meaningful and easier to solve we make the 
following transformations 3 : 



dr d(j) dr J 
d(p dr d(p r 2 



1 du dr du , 1 . 

r ^" = r • s = 3* ' (4 - 28) 



Therefore, Eqn. (|4.26j) becomes: 



* = ( * ) .. ( 4 ) . +(1 _^ )( «£^£ + 4 



Which can be written as: 



(4.27) 



" (- ~ 2 ) 2 ■ {Ju 2 ? + (1 - 2M M )( 24Ca ^ 2 MM3 + JV) 

— ) 2 - J 2 + (l-2Mw)( ^4 + JV) 4.29 

etc? mf 



'^) 2 = 2M 2 Au 4 + (2M-MA)u 3 -u 2 + ^- (4.30) 



where we have defined the dimensionless constant: 



A= «^4 = J^, (4.3!) 



2 We must note that E and -B^ are not equal, E is the energy from the classic relativistic 
orbit equations, while Ek is the quantum energy of the photon 

dr 



3 Using J = r 2 |f 
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In the last form we have used (without proof) the relation E k = V3E 1 which 
will be proven in Sec. 14.2.1 1 Eqn. (|4~68l) . 4 

Vertical Polarization 

Similarly, for the vertically polarized case, using Eqns. (14.181) and (I4.19J) in 
Eqns. (l4~TI)- (l4~ojl to give: 

= r 2 - (1 - — ft 2 + r 2 (l - — )6 2 

2M, „ 8c Q/ 2M,3E?M 
+ (1 ) r ^__|(i (4.32) 



r r- 



t + — (1- — )"Vt (4.33) 



M . 2Af\ -i „ , 2Af 
= f-— (1 ) _1 r 2 -r(l W 

ry £i ry ry 

M . 2M N . 9 4c a 9E 2 M 

+ — i * -Hr^r- 434 



= + — (4.35) 

r 

which, in a similar way as before, gives us the equation of motion for the 
vertically polarized photon: 

du 2 = _ 2M 2 Au A + / 2M + MA)m 3 _ ^2 + ^ / 43 gx 

d<p J z 
4.1.1 Quantum Modified Critical Orbits 

Now the general equation for the quantum modified circular orbits is: 

A 2 = ±AM(2Mu - l)u 3 + 2Mu 3 - u 2 + = /(«) (4.37) 
a<z> .D^ 



4 c Q is dimensionless, while and m e have dimensions of length and inverse-length 
respectively 
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where + is for planar polarization in the r direction, — is for vertical polariza- 
tion in the 9 direction and D is the impact parameter. As A — > Eqn. (|4.37j) 
tends to the classic orbit equation in general relativity, Eqns. (|2.25j) . There- 
fore, in order to determine the magnitude of the quantum correction we can 
calculate the order of A using typical values for m e , D and c a , in Eqn. (|4.31j) . 
Using 5 m e c/h ~ 10 11 for the electron mass (as it is given as inverse length), 
c a = a/3607r ~ 1CT 6 , and the mass of the sun inserted into the critical im- 
pact parameter: D = 3\/3GMq/c 2 ~ 10 2 (given in terms of length), this 
then gives us 6 : 

72r ID" 6 
= JJ^ £U 32 

m 2 e D 2 (10 n 10 2 ) 2 1 ^ 

With the order of A being so small the correction in Eqn. ()4.37|) will be tiny 
compared to the size of the orbit (r = 3GM/c 2 ~ 10 2 ); therefore the modified 
orbits will not differ from the classic critical orbit, given by Eqn. (|2.25l) . by 
very much. 

We will now determine the quantum modified critical circular orbits. In 
order to solve Eqn. (|4.37jl . we can use the fact, from the polarization rule, 
that the modified orbits should be shifted above and below u = -^j by equal 
amounts depending on polarization. So we expect a solution of the type 
u = g|^±(5M; which means we can try a simple modified solution of the form: 

u = Uq + kui (4.39) 

where Uq = ui is the quantum modification, and k is a small constant 

5 As we were working with G = c = h = 1, we must reintroduce these constants to 
obtain the correct order of A 

6 This result is also shown in 
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depending on the quantum modification A. 
Planar Polarized Critical Orbit 



Working with Planar polarization, we can substitute the solution 1)4. 39j) into 
the derivative df/du of Eqn. ()4.37|) : 

-£■ = 2AM 2 u 3 + 3AM(2Mu - l)u 2 + QMu 2 - 2u = (4.40) 
du 

and as A <C 1 and k ~ A, then the only terms of relevance are the ones first 
order in k and A, everything else can be assumed to be ~ 0. Therefore, we 
have 7 : 

2AM 2 u% + QAM 2 u 3 - 3AMu 2 + QMu 2 + \2Mku 1 u^ - 2u - 2ku x = 



AM 9 

Now, substituting this into the trial solution (I4.39J) . we have: 

AM , 

u = u 1 + -^ru ) 4.42 
o 

where u = is the classic solution. Therefore the classic orbit u is 
modified by ^-u to first order in A. Also, with this orbit modification we 
require an associated, modified, impact parameter, which should take the 
form: 

1 1 

'Q 

where 1/D% = 1/27M 2 . The modified impact parameter can be found by 
substituting (|4.42|) and (|4.43jl into (I4.37jl and solving for j3. Doing so, we 



W* = Dl + e (443) 



'u = uq + kui u = Uq + 2ku\Uo u = Uq + ikuiu^ 
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find 



8. 



{^■f = 2M 2 Au i -AMu 3 + 2Mu 3 -u 2 + -^ + (3 = 

do Dn 



2M 2 Au\ - AMul + 2Mu 3 + M 2 Au 



u o 



1 3 1 



" 3 AMn o + ^| + /? (4-44) 

We have 2Mu^ — uf] + ^ = 0, as this forms the classic equation of motion 
for circular orbits. Therefore, 

2M 2 Au\ - AMu\ + ^f 2 ^ - -AMul + P = ® 

3 



Now, as mq = tit, we have: 



3M' 

n AM o 

P = -g-«g (4-45) 

Therefore the modified impact parameter, for planar polarization, is: 

1 1 AM „ 

= If + — "» (446) 

Then, substituting for D , we have: 

1 1 AM 3 _ u 2 Q AM 3 

~ 3(3Mp + ~ T + ~1T U ° 

v 2 

= Ji(l + AMu ) (4.47) 
Vertical Polarized Critical Orbit 

Doing the same for the vertically polarized photon, i.e. by using: 
df 

-=- = -2AM 2 u 3 - 3AM(2Mu - l)u 2 + QMu 2 - 2u = (4.48) 
du 



8 we, again, work to first order in A: u = uq + A^-Uq, u 2 = u'q + A^-Uq, u 3 = Uq + A^-u^ 
and u 4 = ut + A^ful 



CHAPTER 4. QUANTUM MODIFIED TRAJECTORIES 40 
and substituting the trial solution (|4.39jl we find: 

-2AM 2 ul - 6AM 2 u 3 + 3AMu 2 + 6Mu 2 + 12Mk Ul u - 2u - 2k Ul = 



AM 9 



AM > 

u = u (l —u ) (4.50) 

which is equal, but opposite in sign, to 1)4.41)) . as is expected from the po- 
larization rule. Also, as before, the relevant impact parameter is given by 
substituting 1)4.50)1 and (J4.43)) into the negative equation of 1)4.37)) and work- 
ing to first order in A. 

^ = -2M 2 Au 4 + AMu 3 + 2Mu 3 -u 2 + -^ + (3 = 
d<p Dq 

= -2M 2 Auq + AMuq + 2Mu% - MAuq - u 2 

1 , 1 

\ 

Eliminating terms and rearranging, as before, we find: 



3 AMul + ^ + (3 (4.51) 



AM n 

(3 = -—«•» (4.52) 

Therefore the modified impact perimeter, for vertical polarization, is: 

1 1 AM 3 u 2 AM , 



2 1(r>M\2 Q "0 



L> 2 3(3M) 



Wn = — — «r 



?y 2 

^(1-Mmo) (4.53) 



The Modified Orbits 



We now have the circular orbit solutions for Eqn. ()4.37)) and the relevant 
impact parameters: 
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Figure 4.1: Orbit shifts about the classic critical orbit Uo = 1/3M, with 
corresponding shifts in the speed of light (Not to scale) 

• Planar (r) polarization solution (c<l) 

u = u Q + A{^rul) -^ = ^(l + AMu ) (4.54) 

• Vertical^) polarization solution (c>l) 

u = u -A(—ul) - = ^(1-Mn ) (4.55) 

which are displayed in Fig. 14.11 We can note that as the constant A is of the 
order 10~ 32 these modifications are extremely small. 

4.2 Quantum Modified General Geodesies 

Using the solutions for the critical circular orbits and the associated impact 
parameters we can now study how the general photon trajectories are modi- 
fied due to quantum corrections. In the classic case when a photon comes in 
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from infinity, with the critical impact parameter, it tends to the critical orbit, 
as was shown in Fig. 12 .3j and if we slightly decrease the impact parameter 
the photon spirals into the singularity, Fig. 12.51 We will now construct a gen- 
eral quantum modified equation of motion and determine how the geodesies 
change, depending on polarization, as they tend to the critical orbits. 

4.2.1 General Vectors 

From the quantum modification term in Eqn. (|4.1J) we can see that in order to 
construct a general quantum modified geodesic equation we require general 
photon polarization and wave vectors. The wave vector, k^, can no longer 
be represented by a simple constant vector pointing in the direction, as 
given in the orthonormal frame. And, even though the vertical polarization 
will remain a constant, as before, the planar polarization will now be a more 
general vector, constantly changing as the photon moves through the plane. 

General Wave Vectors 

Our general wave vector will be of the form: 

*- = <££.f.£> («•> 

which respects the light cone condition: 

= g^krhP = Ft 2 - F~ x r 2 - r 2 2 (4.57) 

where g^ u is the Schwarzschild metric. Using previous results, of Eqns (I2.12B 
and (|2.13j) . and the fact that we are working in the plane, 9 — |, we can 
represent three of the wave vector components as: 
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where we have defined: 



3 - ^) («0) 



Now using Eqn. 1)4. 5 7|) we can write the final component as: 



% = ^-^=^(1-^,1 (4,0) 



where we have used D 2 = Now the general wave vector can be written 



as: 



Also: 



k^ = (EF-\E(l-^f)^0^) (4.61) 



k, = (^-r^l-^M.-J) (4-62) 



We now have: 



= ^ - ^ 2 (! - ) - ^ = 
as required. Therefore, eliminating -E 1 , we have the general wave vectors: 

* = (1-^.0.^) (4-63) 

*, = ^(l.-F-Hl-^'.O.-D) (4.64) 

We now need to normalize these vectors so, as they come in from infinity 
and tend to the critical circular orbit, they correspond to the critical orbit 
vector (|4.8jl . However, as (|4.8jl is given in the orthonormal basis, and we 
are now working in the coordinate frame, we must use the tetrad (|A.44j) . 
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given in Appendix B, to transform (|4.8j) to its coordinate frame equivalent. 
Therefore, (|4.8j) in the coordinate basis is given as: 



-i\n 



( E k \ 



\E k ) 



E k 



( F~\ \ 



V ~r J 



(4.65) 



This now corresponds to a wave vector for a circular orbit in the coordinate 
frame. To represent the critical orbit we simply substitute for r = 3M, in 
which case F — 1 — — — > |, thus ()4.65|) becomes: 

^ = ^(^3,0,0,^) (4.66) 

Now, if we evaluate ()4.64jl with r = 3M and D = 3a/3M we have: 

F = ^V3(V3, 0,0,^) (4.67) 

Therefore, ()4.67|) is similar to (|4.66j) up to a constant of normalization, given 



as: 



E k = ^V3 = EV3 

Then, the normalized form of (|4.64|) is given by using ()4.68|) : 

, ,, Ei. . „ 1 , D^F ,i D . 

fc" = -4 (F" 1 , 1 — — 5 ,o,- 



General Polarization Vectors 



(4.68) 



(4.69) 



Now we need to construct the, planar and vertical, polarization vectors, a M , 
of the photon; which must be spacelike normalized as: 



-1 a fl k fl = 



(4.70) 
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As before, as we are working in a plane, the vertical polarization, a^, will be 
a constant, and can simply be written as: 

< = (0,0, 1,0) (4.71) 

To normalize this we do as follows: 

ai M = ^ai = -r 3 (0, 0,1,0) ata^ = -r 2 (4.72) 

Therefore, the normalized vertical polarization vector is given as: 

< = (0,0,^,0) 

a lM = -r 2 (0,0,-,0) (4.73) 
r 

These now satisfy both the conditions in (|4.70j) . The planar polarized vector, 
, is given in the plane of r and 0: 

o£ = (0,40,S) (4.74) 

Now, using the two conditions in (|4.70jl we can determine A and B. From: 

^ = 2 M X = (0, -AF~\ 0, -r 2 B) (4.75) 

and the vector ()4.69|) we have: 

k»a 2lx = (1 - ^)I(AF" 1 ) + (^)(r 2 B) = (4.76) 

= A 2 F~ 1 + 5 2 r 2 = -1 (4.77) 

Solving these for A and B: 

F- 1 , D 2 F i „ 
fi = — -(!--— )M (4.78) 
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Therefore, we have: 



A = 2^-1 ™ir- (4-79) 

(1 + ^(1-^)) 1 ' 



„ DF 1 D 2 F i 
A = B = —(l —3 4.80 



and the planar polarization vector becomes: 



a£ = (0,— ,0,-i(l-^)i) (4.81) 



ry /y ' ^"2 



We now have the required polarization vectors: 

1 



r (0,0, 1,0) (4.82) 

ai„ = -r 2 (0,0,-,0) (4.83) 
r 

< = i(0 ) DF,0 ) -(l-^)i) (4.84) 

a 2 » = (0,-^,0,(1-^)4) (4.85) 



where subscript 1 and 2 are vertical and planar polarizations respectively. 
These now satisfy the conditions in (14. 70 J) with the wave vector (|4.69j) . 
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4.2.2 Quantum Modification 

Having derived the general polarization and wave vectors, we can see from 
Eqn. (|3.15|) that the quantum modification given by: 



Sk(a) = -^R abcd k a k c a b a d (4.86) 

also requires the Riemann tensor components in the coordinate frame. In 
Appendix A we have calculated the required components as: 

_2M MF" 1 MF- 1 

n trrt 3 n err6 ^rrc/) 

R^ m = 2Mr Ri em = -M£ R '^ t= _M£ (4.87) 

Using this information we will now determine the form of the general quan- 
tum modification; and from the polarization rule, this quantum correction 
should satisfy the condition: 5k\ = —5k 2 for the two polarizations. 

For vertical (9) polarization we have, by using (|4.69|) and (|4.82j) in ()4.86|) : 

R a bcdk a k c a b l af = R a g c g—k a k c = Rtete^k t k t + R r g r g—k r k r + R^g^g—k^k^ 



MF I F~ l E k F- l E k , 



MF^ ( l Jl-^E k Jl-^E k 
, 1 , , DE k . , DE k . 
M F.ji' ' MElF-\l - ^f) 2ME 2 k D 2 



3r 3 3r 3 3r 5 

ME 2 D 2 



(4. 



Similarly, for planar (r — plane) polarization we have, by using (|4.69ft and 
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m 

-Rabcd^ ^ ^2^2 — -^avcv^ ^ ^2^2 — -^avctp^ ^ ^"2^"2 "^ci^cr^ ^ ^2^2 

2M F~ l E k 2 DF 2 MF- 1 DE k 2 DF 2 



n2 p 

MF" 1 , DE k ..E k L D 2 F„ ,DF S 



7-1 /i D 2 F /i D 2 F 

MF -i F^l- — ng fc yl-— DF, 



r v 7 ^ 



+ 



JJF-'.V-^Vl-'? u 



r v 7 ^ 

'2 



MElD 2 



(4.89) 



Therefore, the quantum modifications, k 2 = 5k (a), for the two polarizations 
di and 02 (vertical and planar respectively) are: 

5k( ai ) = -(%) ME P 2 Vertical (4.90) 
5A; (a 2 ) = (%) M ^ 2jP2 Planar (4.91) 

where — 5k(a\) = 5k(a 2 ), as required. Now, using k^k^ — 5k(a)=0 we can 
write the general equations of motion for the two polarizations: 
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As before, substituting for t, and transforming r — > -, we have: 

We can also write u as a function of time: 

(f ) 2 = (^) 2 (^ - ^ 2 ± (^fV« 5 ) (4-93) 

where we have used the substitution for A given in equation (14.31(1 . Now, 
the Eqns. 14.921 and 14.931 are the general equations of motion, + for vertical 
(6) polarization and — for planar (r-0) polarization. In these equations we 
not only use the impact parameter of the form 1/D 2 but also -D 2 , therefore 
the parameters for the two polarizations are given as: 

1 =% 1±AMu0) ^ D ^^ + O(A) (4.94) 

4.2.3 General Trajectory to the Critical Orbit. 

Now that we have the general orbit equation (14.92(1 we can first test whether, 
for the critical impact perimeter D, the equation ^ — > for first order 
in A. Substituting the critical orbits and the impact parameters given in 
()4.54|) . ()4.55)) and (|4.94jl into ()4.92jl and expanding up to first order in A we 
have: For planar polarization. 

,du^ uZ, A s , AM \ \ r / / AM Nl0 

(^) 2 = -f(l + AMu )-[l-2M(u (l + —u ))[(u (l + —u )] 2 

,ut(l+AMu ) AM ... .... AM AM 5 



■)[1 - 2M[ Uo (l + -—uo)]}[uo(l + -^u )}- 
o o 



(4.95) 
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and similarly for vertical polarization. 

(^) 2 = ^(l-AM Mo )-[l-2MK(l-- r M )]]K(l-- r «o)] 2 

,«JMMua)^ M Ul nu M AM AM 5 

- (— jp )[1 - 2[(txo(1 - -g-«o)]][«o(l - — u o)f 

-> (4.96) 

Expanding these and eliminating all terms of order A 2 and higher, we find 
that the right hand sides go to zero, as required. Thus for the appropriate 
impact parameters these equations behave as they should. 

The next step is to solve equation ()4.92|) for the two polarizations. This 
is most simply done using numerical methods in Mathematica. As a guide, 
we know our solution will be of the form u = u + kui, where u will be 
the classical solution (|2.39j) . ui will be a small modification that pushes the 
critical orbit up or down depending on photon polarization, and k will be 
some constant that is first order in A, i.e. of the form k = As, where s 
will be some number given by the boundary condition: for <fi — > 00 then 
u{4>) — > Uo(l ± ^ uq). If we substitute for the critical impact parameter D 
from (|4.54l) and (|4.55j) depending on the polarization, and then transform to 
u — > uo + Asui, and expand to first order in A, we have non-linear first order 
differential equations in u , u\ and 0. 

• For planar polarization: 

du\{4>) 1 



27M 3 u ((j>) t 



+ 54M 4 n o (0) 6 -6su o (0)n 1 (0) + 18sMu o (0) 2 ni((/))] (4.97) 
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• For vertical polarization: 



dui((j)) 
d(j) 




54M 4 n o (0) 6 - 6SU0OKO) + 18sMwoO)VO)] (4.98) 



These can be solved in one of two ways, (i) is to substitute the classic so- 
lution (|2.39l) for uq and solve analytically, (ii) is to solve and the classic 
equation for ^ simultaneously using numerical techniques. We attempted 
to use method (i) to derive an analytic solution for «i(0), however, due to 
the complex nature of the equation we proceeded to use method (ii), that is 
solving by numerical methods. To do this we set the constant s = 1, and 
then when the numerical values of U\((f>) were determined we could deter- 
mine the constant s so that u(<f>) coincided with the modified circular orbits 
given in (|4.54|) and (|4.oojl . This technique was used for reasons of conve- 
nience, because solving (|4.97j) and (I4.98J) for various values of s would be 
time consuming as each numerical calculation takes a significant amount of 
time; therefore solving them once and then scaling the solution is a more 
convenient method. The results of the equations were plotted as: 



where the constant s was picked to satisfy the condition: — > 00 =>- u(4>) 



polarization. In this way the constant s was determined to be s = 1/3, 
which was tested for various values of M. We have plotted the results of the 
numerical calculation in Figs. 14.21 and 14.31 In Fig. I4.2[ you can see that the 
general critical orbits follow a classic style path, however the orbit splitting 



u(<j>) = u Q (4>) + Asu^cj)) 



(4.99) 




Uq) i.e. the trajectories tend to the critical orbit, depending on 
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Figure 4.2: The quantum modified critical orbits follow classic type paths. 
We have used a factor A = 0.00009 in order to make the quantum corrections 
more visible, where a real value should be of order ~ 10~ 32 , Eqn. 14.381 The 
splitting of the orbits can be clearly seen in the Fig. I4.3| given below. 

is not clearly visible. But in Fig. 14.31 we have plotted a closer view of the 
critical orbits, and here the splitting is highly visible. It can be seen that the 
general trajectories for the planar and vertically polarized photons tend to 
the relevant critical orbits. 
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Figure 4.3: The quantum modified trajectories show clear splitting as they 
tend to the critical points. We used the constants M — g and A = 0.00009. 

4.3 Quantum Modification and the Event Hori- 
zon 

In Chapter. [2] it was shown that when we decreased the impact parameter 
from the critical value (D = 3^/3M to D — 1/10) the photon trajectory, 
given as u(<p), spiraled into the singularity, Fig. 12.51 and when we represent 
the trajectory as a function of coordinate time, u(t), it tends to the event 
horizon, u = 1/2M. From the horizon theorem it was seen that quantum 
modifications have no effect on photon velocities directed normal to the event 
horizon. So this implies when a photon tends to the event horizon at an an- 
gle, e.g. with an impact parameter D = D crit i ca i — 1/10, then the component 
of velocity normal to the horizon should be unchanged, while the component 
parallel to it is modified according to the quantum correction; this modifica- 
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tion should then result in a shift of the photon trajectory, but the horizon 
should remain fixed. In order to test this we used Eqn. 1)4. 92j) with an im- 
pact parameter D = D critica [ — 1/10 to show that the photon trajectories 
still fall into the singularity. We then used Eqn. (|4.93j) . with the new im- 
pact parameter, to study the behavior of the trajectories around the event 
horizon. 

4.3.1 Trajectories to the singularity 

In order to construct quantum modified trajectories, which go past the crit- 
ical orbit and fall into the singularity, we require the impact parameters: 



1 1 V3 y/3AM 1 t . . 

— = F ► D = — — =f — + 0(A) 4.100 

D 2 ( . ^3 _± )2 u 2 10 V ; V ' 

K u oy /(l±AMu ) !0 ; 

where, as before, + is for vertical polarization and — is planar polarization, in 
1/D 2 . For these impact parameters we numerically solved Eqn. ()4.92|) . and in 
Fig. 14.41 we can see that all the trajectories follow a classic type path into the 
singularity 9 . However, near the singularity you can see the splitting of the 
orbits as they tend to u — > oo. In Fig. 14.51 we have shown a magnified view 
of the point where the trajectories cross the event horizon. In this figure 
it can be seen that the planar polarized photon (c < 1) crosses the event 
horizon at a point before the classic trajectory and the vertically polarized 
photon (c > 1) crosses it at a point after the classic trajectory. This makes 

sense, as the planar polarized photons are pushed towards the black hole 

9 The quantum modifications to the classic trajectories are very small, and even if we 
use the hugely exaggerated value of A = 0.00009, as was used in Figs. 14.^1 and 14.21 the 
modification is hardly visible. So, in order to magnify the quantum correction even more 
we used A = 0.0009. 
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Figure 4.4: The quantum modified trajectories follow similar paths, as the 
classic trajectory, to the singularity, with M = ~ and A = 0.0009. 

and vertically polarized trajectories are pushed out, the vertically polarized 
ones must spiral further around the black hole to reach the event horizon 
compared to the classic or the planar polarized trajectories. 

4.3.2 Fixed Event Horizon 

To demonstrate the fact that the event horizon remains fixed at u = 1/2M we 
applied the impact parameters given by (|4.100|) to Eqn. (l4.93j) . Again, solving 
numerically we found that the trajectories tend to the event horizon in the 
classic way, Fig. I4.6| however, they are again slightly shifted. In Fig. 14.71 
a close up of the point where the trajectories tend to the event horizon at 
u = 1/2M, you can clearly see that the quantum modified orbits tend to the 
horizon before the classic orbit. This can be understood by the fact that the 
planar polarized trajectory is pushed towards the black hole, hence it has less 




of a distance to propagate before it reaches the horizon, and even though the 
vertically polarized trajectory is pushed outwards, the fact that it has a faster 
velocity than c = 1 it reaches the event horizon before the classic trajectory. 
The other more important thing we can note from this calculation is that the 
event horizon remains fixed at the classic point u = 1/2M, due to the fact 
that, as was stated earlier, the quantum correction has no effect on photon 
momentum normal to the horizon. This means even though the quantum 
correction implies velocities greater than the speed of light, no photons can 
escape from the event horizon, so the black hole remains black. 
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Figure 4.7: The classic and quantum modified trajectories reach the event 
horizon at differing times, with M — | and A = 0.0009. 



Chapter 5 

Quantum Modified Schwarzschild 
Metric 

Now, using our previous results, we can construct a new metric that en- 
compasses the quantum corrections due to vacuum polarization. This metric 
should then be a sum of the classic Schwarzschild metric and the polarization 
dependent quantum corrections we derived in Eqns. 1)4. 88|) and (|4.89|) : 

= (<j F ~^aV)Fr = (5.1) 

5.1 Construction of the Metric 

Using the results in Sec. 14.2.21 we can write the first components of the quan- 
tum modified metric Q\ t and Qf t , for vertical and planar polarization respec- 
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tively, as: 



9u 
F - 

F - 

9tt 

F - 
F - 



a P? n h n d 
m e 

*?<-^)<±) 

mi r r z 



AMD* 



9 



-Fu 3 



-RtbtdO^a^ 



8c a 2M DF 2 



,MF 



AMD 2 



)( 



9 



(3D 2 F 2 u 5 - Fu' 



D 2 F 



(5.2) 



(5.3) 



where, in the last lines of each component, we have used Eqn (|4.31j) to replace 
the constants with A, and we've made the transformation r — > l/u. Now, 
doing this for all the other components we can construct the metrics for 
vertical and planar polarizations: 



Vertical Polarization: = u(0, 0, 1, 0) 



( F — BFu 2 



\ 
















J_ i Bu 3 
F F 

o -4 





Planar Polarization: a% = u(0, DF, 0, L>g) 
/ F + BFu 3 (l -3D 2 Fu 2 ) 









-\-Bu J 



Q 



fill 
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F F 
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(5.4) 
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where B = These are now the relevant quantum modified Schwarzschild 

metrics for vertical and planar polarizations 1 . 



5.2 Dynamics with the Quantum Modified Met- 
ric 



We can now derive the same equations of motion as ()4.92|) . but now we can 
do it simply with the quantum modified metrics. Using the general wave 
vector in a plane: 

k" = (t,f,<i,0) 

= ^'■-4 M " 2) (5 - 6) 

where we have transformed to r — > 1/u, and we substitute for i and (f) as 
before. Applying this wave vector to (|5.4j) and (|5.5l) we find: 



= Q^W 

- 2= ^- Fu2+( ^ )Fu5 (5 - 7) 

for vertical polarization, and 



- O 2 = h - f - 2 - <^« 5 < 5 ' 8 > 

for planar polarization. These are identical to the equations we derived in 
Sec.0~2H 



1 In both the metric, and vector, 0,%, we have made the substitution r — » 1/u and 
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Radial Geodesies 

We can now show that the metric for vertical polarization is consistent with 
the fact that radially projected photon trajectories are not modified. By 
using the vertical polarization metric, 1)5. 4j) . and a general radial wave vector, 
(&" = (t,f,0,0), we find: 



o = 

dt 1 , . 



which is identical to the classic radial geodesic equation, 1)2.30)) . However, 
the same is not true for the planar polarization metric. Due to our deriva- 
tion of the polarization vectors in Sec. 14.2.11 we constructed a very general 
vertical polarization vector and then normalized it; however, the one for pla- 
nar polarization was constructed for the case where ^ 0, as can be seen 
in Eqns. 1)4.60)1 and 1)4.76)1 : therefore the planar polarization vector has 
dependence "mixed" into it through the substitution of <j>: 

.dr . 2 2 .dt. 2 2 .d(f>. 2 



(5.10) 



We can note that it's a result of the substitution ^ = fj, in the polarization 
vectors, that we acquire an extra parameter D in our quantum correction 
(apart from the one introduced through normalization, which was incorpo- 
rated into B = AD 2 / 9). We can then say that if we need to use the metric 
for a radial trajectory with ^ = we can just set D = 0, as this would 
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lead to the removal of the dependence. In this way, the planar polarized 
quantum modified metric for radial trajectories is: 



( F + BFu 3 






\ 

F F 

o -4 o 



(5.11) 





-\ o 

/ 

Then applying this to a general radial wavevector, we again find the classic 
radial geodesic equation: 



= G(<P=0)l u k»k» 
dt 1 
dr F 



(5.12) 



Therefore, we have metrics for the Schwarzschild space time that incorpo- 
rate quantum corrections due to vacuum polarization; and these metrics are 
consistent with classic results. 



Chapter 6 
Summary 



For Schwarzschild Spacetime we showed, in the orthonormal frame, that, 
due to quantum corrections, the stable circular orbit at uo = 1/3M splits 
depending on the polarization of the photon, and the new modified circular 
orbits are given by the classical orbit plus a correction term to first order in a 
dimensionless constant A: u = uq+Au\, where A = 2 j^ fc and has an order of 
10~ 32 . As stated by the polarization sum rule this splitting of the critical orbit 
is equal in magnitude but opposite in sign for the two polarizations. We found 
that the vertically polarized photon (c>l) is pushed out by a correction of 
«i = — ^Uq, while the planar polarized photon is pulled in by a correction of 
u i = ^ u 1- It was a l so found that this orbit shift also requires an appropriate 

2 

modification of the classic impact parameter -kz = ^r. This modification, 
for the quantum corrected orbits, took the form 752 = + ADi, and again 
as for the orbit shift, the change was equal in magnitude but opposite in 
sign for the two polarizations: D\ = for planar polarization and D\ = 
_Muq vertical polarization. Using this information, for the splitting of 
circular orbits, we then constructed the quantum corrected general equations 
of motion for the Schwarzschild spacetime. Using these equation of motion 
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it was shown that a photon starting at u = 0, with the appropriate critical 
impact parameter, tends to the critical orbit associated with that impact 
parameter - and the trajectory follows a similar path to the classic case. 

We then went on to show that, using the general quantum corrected equa- 
tions of motion, a photon projected towards the black hole with an impact 
parameter less than the critical value falls into the the singularity, in terms 
of the angular distance (0) . Although the photons follow a classic type path 
into the singularity, the trajectories are slightly shifted according to polar- 
ization. The planar polarized photon crosses the horizon before the classic 
trajectory, and the vertically polarized one crosses it after, this corresponds 
to the fact that planar polarized photons are pulled towards the black hole 
and vertically polarized ones are pushed away, hence the vertically polarized 
ones need to go a further angular distance to reach the event horizon. 

In terms of coordinate time (t) we found that the photon trajectories 
tend to the event horizon. Therefore, the quantum corrections do not shift 
the classic event horizon from u = ^u, which corresponds to the horizon 
theorem. Also, we found that, although the quantum corrected orbits follow 
a classic type path to the event horizon, the point at which they hit the 
horizon is, again, slightly shifted depending on polarization. However, this 
time both polarizations hit the horizon before the classic trajectory. The 
planar polarized photon tends to arrive at the horizon first , then the vertically 
polarized one, and finally the classic photon. This could correspond to the 
fact that the planar polarized photon, although it has a velocity lower than 
the speed of light, has less of a distance to go, as its trajectory is pulled 
towards the black hole. For the vertically polarized case, even though it has 
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a faster than light velocity, it has a further distance to go to reach the horizon 
as its trajectory is pushed away from the black hole. 

Having determined the equations of motion, with the quantum correction, 
we then went on to construct a Schwarzschild metric that incorporates the 
quantum correction: Q^ v = (g^ + g^), where the correction: g^ v was again 
first order in A. We showed that with this metric and a general photon 
wave vector, we obtain the quantum modified equations of motion, as 
before. Also, this new metric confirms the horizon theorem, that is, when 
we use a wave vector indicating a radially projected photon we obtain the 
classic equation of motions. This was fine for vertical polarization, however, 
in the planar polarization case we had a problem; we had previously used 
a substitution that mixed a "hidden" radial angle, 0, into our polarization 
vector (as in general orbits the planar polarization depends on (ft). By tracing 
back to the origin of this substitution we found that, in order to study radial 
trajectories, we need to set D — 0, this then removes the (ft dependency; 
this then gives us the classic equation of motion for planar polarized radially 
projected photons. 

So in conclusion, after studying the dynamics of null trajectories in Schwarzschild 
spacetime we derived the polarization dependent photon trajectories to first 
order in the constant A (which is dependent on the fine structure constant, 
the mass of the star, mass of the electron, and the energy of the photon). We 
then incorporated these modification into a general quantum modified metric, 
which could also be used to derive the general quantum modified equations 
of motion. Also, the results of this work coincide with the conditions of the 
horizon theorem and the polarization sum rule. 



Appendix A 

Schwarzschild Spacetime 



The 'Schwarzschild geometry' describes a spherically symmetric spacetime 
outside a star, and its properties are determined by one parameter, the mass 
M. The Schwarzschild metric, in spherical polar coordinates, takes the form: 



9w 
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— r 2 sin(#) J 



(A.l) 



This type of spacetime geometry is said to be static, due to the fact that (i) all 
metric components are independent of t, and (ii) the geometry is unchanged 
by time reversal, t — > —t 1 . 

A.l Geodesic Equations 

The equations of motion in Schwarzschild spacetime can be derived by using 
the covariant form of Newton's second law of motion, 

dx v 



Dt 

dr 



_|_ n p 



vp 



dr 



P 



(A.2) 



X A space time with property (i) but not necessarily (ii) is said to be stationary i.e. a 
rotating star/black hole 
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where we have p M = ?ri%-, and for free fall we have F M = 0, then the 



equations of motion are: 



d 2 x fl dx v dx p , , , 



Now, using this and the appropriate Christoffels we are able to determine the 
equations of motion. Alternatively, by using the Schwarzschild Lagrangian 

L = ~[(1 - — )dt - (1 - —)- x dr 2 - r 2 d0 2 - (r 2 sin 2 fl)c/0 2 ] (A.4) 



2' 



we find 
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(A.5) 



(A.6) 



(A.7) 



(A.8) 
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and substituting into the Euler equation 

dL d dL 



dx^ dX dx^ 

we have the equations of motion for Schwarzschild spacetime 



(A.9) 



2M S - 2M . , A , 

= (1-— )t+— ft = (A.10) 



= — t 2 -r0 2 -r0 2 sin 2 + (1 ) -1 r 

= -r 2 sin0cos0</> 2 + 2rr0 + r 2 (A.12) 
= 2rr sin 2 00 + 2r 2 sin 6 1 cos 6<j> + r 2 sin 2 9(f> (A. 13) 



A. 2 Riemann Components in the Orthonormal 
Frame 



The spacetime line interval, Eqn. Ijl.llj) . can be rewritten using tetrad trans- 
formations (discussed in the next section): 

u t = (l-—)^dt u r = (l-—Y^dr u e = rd9 to* = rfan(9)d<f> 

(A.14) 



ds 2 = (J? - U) 2 - (u 9 ) 2 - (u^) 2 



(A.15) 



This now has the Minkowski metric 



9nu 



( 1 \ 
0-100 
0-10 
\ -1 ) 



(A.16) 
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Noting that the metric can be written as: 

dg^u = c*V + w vft ( A -17) 

and using the fact that: 

d9,u = 9 -§^dx« = (A.18) 

we then have|2]: 

U)^ = -LO^ (A. 19) 

which implies to is antisymmetric, i.e. it has only six unique components and 
i^ixv = for \i = v\ and we have the conditions: 

w? = 4 = -<W (A.20) 

We can now write the exterior derivatives 2 of u; M 



1, 2.W, ,1X1 
M, 1M i 



dJ = -(1 — )-2(_ )d r Adt 



r 



2 



1 )~*u; r Au* (A.21) 



dw r = (A.22) 

du e = dr Ad6 = -(1 - — )V A J (A.23) 
gL/ = sin(0)dr A d(p + r cos((p)d9 A d(p 

= V-)^A^ + ^A^ (A.24) 

ry ry ry ' 



Now, using Cartan's equation 

d^ 1 = u a A u£ + IF (A.25) 

2 The exterior derivative: d = dx^-^jz acting on a 1-form A v = A(x)dx u gives dA v = 
dx^ A dx"^ 1 , where dx^ A rfx" = for /i = i/ 
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with zero torsion = 0), we can write 

dtJ = uj t Atul + lu 9 AluI + w+ Alu^ (A.26) 

duf = J- Au r t +u e Aul + L}* Aul (A.27) 

du e = J Auo e t +uo r A uj e r +uo* Aul (A.28) 

du* = uj 1 A uf + uj t A uf + u e A u)q (A.29) 



Comparing Eqns. (TOTTl - lTOl with Eqns. (|A.26D - (|A.29I) we find the six 
unique components of u% as 



, M . 2M, i , M , 

u* r = -rl )-2j = dt = u r t A.30 

J r = -1 )*uu e = 1 )*de = -uj r e A.31 

ty <y qy 

uf = -(1-— )3w* = Bin(0)(l-— )3d0=-wJ (A.32) 

7™ 1™ ^" 



cot(0) ^ = cos(W = (A33) 



r 



^ = ^ = (A.34) 
= oA = (A.35) 



Now, using: 

^ = d< + < A < (A.36) 
we can write the six independent components of the Riemann tensor: 

R l r = dJ r + iJ a A u° = duj l r + J e A uj e r +uj\A cuf 
2M . 2M 



dr Adt = '-^-cu r A uj 1 

2M 



^rrt Rtrrt 



(A.37) 
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Rl 



dJ + (J A = da/ + A + a/ A uol 



n rr g 



2M. i2M , M e 
)"2 ^^rfr Add — —u T A w fl 



— Rf)r, 



M 



6rr0 



(A.38) 



i# = duf + utAuj^duf + uf Aul + uf Au e r 

2M i M 2M i 
= cos(0)(l )3d0 A d<f> + sin(6») — (1 )~5dr a 

2M 

+ cos(0)(l ) 1 2d(PAd9 

r 

cot(0). 2M.i e . M . 
= ^(1 )*uj e A^ + —uj r Auj* 



ry £t ry ry 

M 

3 



cot^)^ 2M 
r 



r 



(1 — ) 5cJ A uj = —u r A u q 

M 



(A.39) 



Rl 



duot + out A u% = duf + out A ui + oof A 



-sin(0)d0 Ad0-sin(0)(l 



2M 



)d(f> A dO 



1 , . 1 ,„ 2M, fl . 2M fl . 
w A — -(1 )uj e Auj^ = -J A 



iy»2 ry 2 " iy» 



-R 



2M 



(A.40) 



du£ + uj a A cu? = dJ t + J r A uj T t + uj , A u)\ 



M . 

f? e - 



2M, i 



|2d6» A dt 



M 



-uj p A a; 1 



r 

—Retet 



M 



(A.41) 
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Rf 



du;? + uoi A u« = dujf + ojf A 0J r t + ujg A J{ 
^(1 - — )a s in(0)cZ0 A dt = ^ru* A J 



R 



t(i>t 



-R 



M 



<pt(j>t 



(A.42) 



A. 3 Tetrad Transformation 



In Eqn. (|A.14|) we transformed from the coordinate frame to an orthonormal 
frame. This is achieved by using the tetrad of the form: 

\ 



/ 

F~$ 

r 

\ rsin(fl) / 



(A.43) 



This takes a 4-vector, given in the coordinate frame, and maps it to the 
equivalent in the orthonormal frame. We can then determine the inverse 
tetrad, i.e. a tetrad that takes vectors in the orthonormal frame and maps 
them to the coordinate frame. We find this by using the inverse of (|A.43J) . 



( F- 



J a 










Fs 
i 

r 











\ 



(A.44) 



m ) 



which has the property: 



(A.45) 



Therefore, (|A.44ft takes vectors in the orthonormal frame and maps them to 
the coordinate frame. 
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A. 4 Riemann Components in the Coordinate 
Frame 

In Sec. IA.2I the unique components of the Riemann tensor were calculated 
in the orthonormal frame. In order to determine the components in the 
coordinate basis we will use the transformation tetrad given in Eqn. rrm 
Before we can us the tetrad, (jA.44j) . we need to adjust it so it is able to take 
covectors in the orthonormal frame, rather than vectors: 

/ Fl \ 



F~* 
r 
\ Or/ 



(A.46) 



which is just Eqn. I) A. 43)1 . Then using this we have: 

=>- ^(coordinate frame) = (e _1 )^y4 a (orthonormal frame) 

Now using this tetrad we can transform the orthonormal Riemann tensor 
components as: 



Rfivap = (e x )^(e l ) b u (e l ) c a (e l ) d pR a bcd 
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Now using the components given in Sec. IA.2I we have: 

Krt = (e- 1 )^(e- 1 );(e- 1 );(e- 1 )X rt = ^^-^-^n-^) = -^ 
R'erre = {e^tie^Xie-X^feRerre = rF-^r(~) = 

Km = (e" 1 )^(e^)*(e" 1 )^(e" 1 )^ = r^rFn-^) = -^ 
Kt,t = (e- 1 )J(e- 1 )J(e- 1 )J(e- 1 )^t = r^r^(-^) = -^ 

(A.47) 

Therefore, the Riemann tensor components in the coordinate frame are given 
as: 

P' - f?' MF ~ 1 Tf MF ~ X 

U trrt — r 3 n 9rr8 ~ f n <Prr4> ~ ^ 

MF MF 

R^aa^ — 2Mr Rotot — ~ R<t>t<t>t ~ ~ 



(A.48) 
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